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We consider an interacting Lifshitz field with z — 3 in a curved spacetime. We analyze the 
renormalizability of the theory for interactions of the form with arbitrary even n. We compute 
the running of the coupling constants both in the ultraviolet and infrared regimes. We show that 
the Lorentz- violating terms generate couplings to the spacetime metric that are not invariant under 
general coordinate transformations. These couplings are not suppressed by the scale of Lorentz 
violation and therefore survive at low energies. We point out that in these theories, unless the 
effective mass of the field is many orders of magnitude below the scale of Lorentz violation, the 
coupling to the four-dimensional Ricci scalar i}"^^ Rcji^ does not receive large quantum corrections 
^ 3> 1. We argue that quantum corrections involving spatial derivatives of the lapse function (which 
appear naturally in the so-called healthy extension of the Hofava-Lifshitz theory of gravity) are not 
generated unless they are already present in the bare Lagrangian. 

I. INTRODUCTION 

As there is no hint about how the laws of physics are modified at extremely high energies, it is reasonable to explore 
the possibility of high-energy violations of the low-energy symmetries and test how robust the predictions in the 
infrared are when departures from those symmetries in the ultraviolet are considered. In particular, in the last years, 
there has been a growing interest in field theories where Lorentz invariance is explicitly broken at high energies [l| . 

There is an additional motivation to consider these kinds of theories, which comes from gravity. In the theory 
proposed by Hoi^ava , one introduces a preferred foliation of spacetime and allows for higher spatial derivatives in 
the Lagrangian. The invariance under general changes of coordinates is lost, but the higher spatial derivatives improve 
the quantum version of the theory. It becomes a power-counting renormalizable theory of gravity if one includes 2z 
spatial derivatives, with z > 3. At low energies, one expects that there will be a region in the parameter space of the 
theory in which the well-tested predictions of general relativity are recovered. 

When matter fields are coupled to Hofava's gravity, in principle, one has two options for their corresponding bare 
Lagrangians: one could consider that they also have higher spatial derivatives, breaking explicitly Lorentz invariance 
in the matter sector, or that they do not have them. The latter was considered in Ref. [1] with the main goal 
of making the model phenomenologically viable in the sense that the observational constraints on Lorentz- violating 
effects (which in this case, are induced only by the gravitational coupUngs) could be naturally satisfied. However, as 
was shown in Ref. (3] , the renormalizability of the theory implies that one must also include higher derivatives in the 
matter sector. Indeed, even for free quantum matter fields in a classical gravitational background, z > 3 is needed in 
order to avoid the existence of divergences with more than two time-derivatives of the metric Q • Scalar fields with 
such higher spatial derivatives, which obviously break Lorentz invariance, are named Lifshitz fields, and are the main 
subject of the present paper. 

Various investigations on the properties and applications of Lifshitz fields with z = 3 can be found in the literature; 
see for instance Refs. [B-Q] in the cosmological context, and Refs. [1, Q for nonperturbative studies in fiat space. 
There is an important question that deserves to be addressed, which is whether one can recover the usual Lorentz- 
invariant field theory predictions at low energies or not, for Lifshitz fields. On this subject, in Refs. [l^, [lli, the 
authors considered models with z = 2 scalar fields (and also models with fermions as a Yukawa like model) in higher 
spatial dimensions with the main goal of analyzing the renormalization group (RG) evolution of the "speed-of-light" 
couplings Ci's for one and more than one field (see also Refs. [l3i[3)- The one-loop renormalized effective potential 
for a Lifshitz scalar field with z = 2 at finite temperature has been analyzed in Ref. [IJ] . One motivation for the choice 
of z = 2 and higher dimensions is to obtain nontrivial results without having to compute high-loop contributions or 
diagrams with many legs. The results in Refs. [l^, El indicate that the recovery of the usual Lorentz-invariant 
theories at low energies is nontrivial, requiring in general a fine-tunning of the parameters in the ultraviolet. The 
requirement of such fine-tunnings seems to be generic in the sense that even thought the different q's can run to the 
same value in the IR, in simple models, the running is in general too slow to satisfy the phenomenological constraints. 
However, in Ref. [isj . the authors have pointed out that a faster running can be achieved in some scenarios containing 
a large number of hidden fields beyond the standard model ones. Clearly, the studies carried out so far are still not 
conclusive, and this question should be further addressed. Although this latter issue goes beyond the scope of the 
present paper, one of our goals is to move forward in this direction by extending nonperturbative RG techniques for 
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z = 3 Lifshitz fields. 

The analysis of quantum field theory with z 7^ 1 is therefore of interest for several reasons. On the one hand, if 
one considers Hofava's gravity as a useful arena to understand different aspects of quantum gravity, it is unavoidable 
to include these kinds of matter fields. On the other hand, this analysis can be considered as a simplified situation 
for the study of the renormalizability of Hofava's gravity beyond power-counting analysis. Finally, one could address 
interesting questions as, for instance, the emergence of the usual Lorentz-invariant theories at low energies. 

With these motivations in mind, in this paper, we study the theory of self-interacting quantum Lifshitz fields with 
z = 3 in curved backgrounds. We will consider interactions of the form with arbitrary n, that are allowed by 
power-counting. Looking at the Heisenberg equations for the quantum field operator, we will discuss the one-loop 
renormalizability of the theory. This turns out to be a rather simple task, since the divergences are those of flat 
spacetime. It is also possible to compute the running of the coupling constants in the ultraviolet regime. In order to 
analyze the infrared behavior of the coupling constants, we generalize a method based on the evolution equation for 
the effective potential to the case of Lifshitz fields. 

We consider the quantum fields in a curved background so as to analyze the couplings to the spacetime metric 
induced by the self- interaction. It is well-known that, in theories with z = 1, it is necessary to include couplings to 
the curvature of the form £,R(jP in order to absorb the divergences of the theory. We will sec that, although these 
kinds of terms are not necessary to renormalize for z = 3, the self-interaction induces finite couplings to both the 
3-curvature and the extrinsic curvature. Somewhat surprisingly, these couplings are not suppressed by the scale of 
Lorentz violation, leaving a footprint of the ultraviolet behavior of the theory in the infrared. We will also analyze 
whether quantum corrections induce terms containing derivatives of the lapse function or not, as needed in the 
so-called healthy extension of Hofava's theory [l^ . 



II. INTERACTING LIFSHITZ FIELDS IN CURVED BACKGROUNDS 

As already mentioned, there are several reasons for considering a z = 3 scalar field in a curved background, amongst 
which is to move forward in the understanding of quantum effects in the framework of the Hofava-Lifshitz theory. 
As was shown in Ref. [l9| , a power-counting analysis similar to the usual one can be applied for any z, provided 
one substitutes the standard scaling dimensions of the operators by their "weighted scaling dimensions" , i.e. by the 
dimensions implied by the assignment [x]w = —1 and [t]w = —z. A simple power- counting argument along the lines 
of Refs. [lO, [2l| shows that, in this case, the interacting theory is renormalizable for an interaction of the form (/>", 
with n arbitrary. 

To review this, let us suppose the action of the scalar field in D -|- 1 fiat dimensions S — Sfiee + S-mt, with a free 
part of the form 



- (/. [m2 + A + • • • + A2-2^ (A)"] 0} dt d^x, (1) 



where A is a parameter with momentum dimensions that controls the Lorentz- violation scale. For the interaction part 
of the action, we consider a polynomial in of degree Umax'- 



p{(j)) dt d^x^ ly -^0" ) dt d^x, (2) 




where each coupling constant A„ has dimensions 

[A„]=il/^+i-«^. (3) 
For a generic Feynman diagram with L loops and / internal lines, the superficial degree of divergence [20| is 

5 = L{D + z) -2Iz = L{D - z)-2z{I - L). (4) 

As / > -L in general, for z = 3, we can expect at most a logarithmic divergence {5 = 0). Setting z = 1 reduces this 
expression to the standard result for Lorentz-invariant theories. Considering, in particular, a given diagram with only 
one type of vertex with n legs, and using Euler's theorem for graphs along with other identities regarding the number 
of 71-legged vertices and the number of external lines E, the expression ^ can be rewritten as 



5^ 



(i5 + z)(l-^)-2. 



l-^\{D + z). (5) 
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Now, for a given process, the number of external lines E is constant, while the number of loops L and internal lines 
/ increase with every order in perturbation theory. Therefore, d will not increase for further corrections if 

{D + z)fl-^j-2z<0, (6) 



which translates to 



n < 



2{D + z) 



D-z 



n < oo, 



z < D, 
z> D. 



(7) 
(8) 



This shows that for z > D, the theory is power-counting renormalizable for any power of (/) in the interaction 



Lagrangian. In the Lorentz-invariant case {z = 1), the highest power allowed is n^naxiz = 1) 



the standard result rimax = 4 in 3 + 1 dimensions. 



2{D+1) 
D-1 



which gives 



We will consider a Lorentz noninvariant action in a curved spacetime. Using the Arnowitt-Deser-Misner formalism, 
the spacetime interval is written as 



^N^dt'^ + g,j{dx' + N'dt){dx^ + N^dt), 



(9) 



where, as usual, Latin indices refer to the spatial coordinates, i,j = 1, 2, 3, N and Ni are the lapse and shift functions, 
respectively, and gij is the spatial metric. The explicit form of the action is given by 



S^= dtd^xNyf^ 



1 



m"0^ - bid,(l)d'(j) 



A2 



(A0)^ 



A4 



E 

p=2 



A 



2p 



(10) 



where the operator A is the spatial Laplacian associated to gij . The sources of Lorentz violation are the coefficients 
bi, ^2 and 63, while, as before, A is a parameter with dimensions of mass that sets the scale at which this violation 
occurs. Furthermore, we included interactions of the type (jP^ up to an arbitrary value of Pmax = nmax/'^, which 
are not renormalizable in the Lorentz-invariant case when n^ax > 4. Neither the kinetic terms nor the interaction 
terms considered here are the most general; for instance, as was emphasized in Ref. [isj . the foliation-preserving- 
diffeomorphism symmetry of the theory also allows the inclusion of spatial derivatives of the logarithm of the lapse 
function N. One could also consider interaction terms containing derivatives of (j). We restrict for the moment to this 
particular choice, and postpone the discussion about more general Lagrangians until Sec. IVI Bl 
The classical equation of motion derived from the action ^TU\\ is 



Ny/g 

where we have defined the operator 



E 

p=2 



2p 



(2p-l)! 



(11) 



b3 



-2^V^A (NA...) + {NA'J + f^^gA^ (NA...) 



2bld, {Ny^dl 



(12) 
(13) 



Proceeding along the same lines as in Ref. [2^, we write the quantum field (j) as its mean value ((/>) plus quantum 
fiuctuations (h around that value, which is 



(</>) 



(14) 



For this decomposition to be consistent, the expectation value of the fiuctuations must vanish {(/)) = 0. Replacing 
Eq. (|14p in the equation of motion (|11[) and taking its mean value {■■■), we are left with the evolution equation for 

^0 = (0): 



(□4 



X - m (po 



Pmax 2p-l 

EE 

p=2 fc=0 



A2p 



k\{2p-k- 1)! 



0, 



(15) 
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while, on the other hand, subtracting this last equation to the equation of motion (jlip . there remains the equation 
for the fluctuation 6: 



NV9 



Pmax 2p-l 



A 



2p 



p=2 k=Q 



fc!(2p- fc - 1 



= 0. 



(16) 



In the 1-loop approximation, the contributions of the expectation values {4>^), (0'^), etc., are much smaller than {4''^). 
Hence, only the terms with fc = and fc = 2 of the sum in Eq. and the term with fc = 1 of the sum in Eq. 
are kept. Then the most important quantum correction to the evolution equation for 00 is proportional to (0^), 



NV9 



p=2 



A 



\2p 



2p+2 iV2\ 



(2p-l)! 



while the fluctuation cj) satisfles a free field equation with variable mass 

iVVff [Ul{x)- M'^ix)] (^(x)~0, 

where 

A^ 



M2(x) = + E 



;\2p ,2p-2 



p=2 



(2p-2)! 



(17) 



(18) 



(19) 



In Eq. (|17p . the sum has been rearranged to group the terms proportional to the same power of 

III. ADIABATIC EXPANSION OF THE QUANTUM CORRECTIONS 

To study the rcnormalization of Eq. (jl7p . it is sufficient to calculate the quantum corrections in an adiabatic 
expansion, i.e. derivatives of the metric. Each order in this expansion is less divergent than the previous one, so the 
full ultraviolet behavior of the corrections is exhibited in the first orders of this expansion. In the case of z = 3, as 
we already mentioned, the only divergence we can get is logarithmic, so in principle, only the zeroth adiabatic order 
correction would be divergent so that the rcnormalization procedure could take place exactly as in flat spacetime. 
Nevertheless, it is interesting to calculate higher adiabatic order corrections; on the one hand, to find the couplings 
with the curvature generated by the quantum effects, and on the other hand, to see whether the usual z ~ 1 case is 
recovered when the appropriate limit is taken. In this last case, it is necessary to go up to second adiabatic order 
when performing the rcnormalization procedure at the level of the evolution equation [22[ . 

The expectation value of (jy^ taken in the vacuum state |0) can be calculated from the Fcynman propagator Gf{x, x'), 



m\xm 



lim Ini G-p{x, x'), 



for which it is necessary to solve the equation for the Green function of the fluctuation, which is 

N^[Ul{x)- M'^ix)] Gf{x,x') =d{x-x'). 



(20) 



(21) 



Throughout this paper, the mass M will be taken as constant. An improvement on this approximation would be 
to consider a derivative expansion on the mass, much like the adiabatic expansion in derivatives of the metric, which 
in this case, would be equivalent to a derivative expansion on the mean value of the field 0o. However, the terms 
including derivatives of (f>o will not be relevant to analyze the rcnormalization. 

We will consider, in turn, two different approximations to simplify the calculation of (0^) up to second adiabatic 
order. First (in this section) is a weak gravitational field approximation, and after that, we will choose an ansatz 
for the metric and perform a derivative expansion (see Sec. IVI A[) . In both cases, the goal is to solve for the Green 
function in the adiabatic expansion. 

We first consider small perturbations around fiat spacetime: 



N = l + Sn, 



9tj 



(22) 



5 



keeping terms up to linear order in these perturbations. The Feynman propagator at zero order in the metric 
perturbations reads 



Gp \x, x') — — 



k-{x — x') 



(2^)4 {u;l-P,-^ey 



where 



c.,. = \/M2 + 26i|fcp+2^|fc|4 + 2^|fc|6. 



The first order contribution, following the steps from Ref. [1], can be written as follows: 



G^p\x, X ) 



,4 „ d^ki d'^k2 d'^p 
d X 



^ik-i-{x — x") pip-x" pik2-{x" —x') 



(27r)4 (27r)4 (2^)4 (^2^ _ k,l) [uj^ - k^l) 



fk2 (P), 



(23) 



(24) 



(25) 



where fkip) is a function of fco j i Po i a-nd pi which is linear in the metric perturbations and can be found in the 
Appendix. To obtain the different orders of the adiabatic expansion, the integrand of Eq. (j25p must be expanded in 
powers of pi and pq. The coincidence limit can be taken trivially and, after Wick rotating to k^ — —iko, the integral 
in ki can be performed easily. After a straightforward but tedious calculation, the following expression for the zeroth 

adiabatic order of the coincidence limit of Gp\x,x') is reached: 



\G 



(l)iAD(O) 



d^k 

(2^ 



Dm 



Sn 



Dl{k) 



DUk)h,,k^k^^ 



(26) 



where 



denotes the coincidence limit, and 

(2)°+l 



Dtik) 



{2n) 



a+l-26 
■'fe 



dx 



-oo (1 



X 



2\b 



l(^ir!Z , ,a+i-2b r(^)r(-^ 

(271-) ^k T(b) 

0, 



a even, 
a odd. 



(27) 



The corresponding expressions for the first and second adiabatic orders can be found in the Appendix. 

Finally, after a few more steps regarding the angular integrals, taking the imaginary part and putting everything 
together, the results for the adiabatic orders of {(j)^) give 



l2xAD(2) _ 



8^2' 

0, 

1 

96^ 



II) {d^h - 2dtd.,5N') + Iq {didjh'^ - did'h) 
- I 2/o 



(28) 
(29) 



(30) 



where we have defined the following integrals: 



Iq 



du ■ 



lo 



du 



u 2 

du 



U 2 d LOf, 



du ■ 



du^ 



'r + l 



du ■ 



d\ 



du^ 



(31) 



and u = k^. By simply power counting, one can verify that, provided 63 7^ 0, Iq is logarithmically divergent, while the 
remaining integrals are finite. Therefore, as expected from power counting, only the zeroth adiabatic order correction 
is divergent (which is the one involving Iq), and hence is the only one we need for analyzing the renormalization 
process (see the next section). We will come back to the result for the second adiabatic order in Scc lVIl when we 
study the quantum corrections to the couplings with the curvature. 
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IV. RENORMALIZATION AND RUNNING OF THE INTERACTION CONSTANTS 



111 the previous section, we obtained the 1-loop quantum corrections to the evolution equation of {(jy) up to second 
adiabatic order. We found that only the zeroth adiabatic order correction is divergent. Hence, the renormalization 
procedure only requires counterterms independent of the curvature, which will cause the renormalization of the 
interaction constants \2p and the field's mass m. 

At high energies, the behavior of the running couplings can be obtained from the minimal subtraction renormal- 
ization scheme. The divergent integral /q is regularized to give a finite expression using dimensional regularization, 
and the counterterms are choosen to precisely cancel the poles in Z) = 3. This process requires the introduction of an 
arbitrary mass scale /i in order to keep the physical couplings with the correct dimensions. 



^2p 



(32) 



where the superscripts B and R stand for bare and rcnormalizcd, respectively, SX2p are the counterterms and e = 3 — D 



2p 

Replacing the bare quantities for the rcnormalizcd ones in the evolution equation, it reads 

p 



2 ^4 



E 



(p-l)e 




'^A2p + 



^2p+2 




{2p-l)\ 



= 0. 



(33) 



The evaluation of Iq is more complex than for the usual case, so we sketch here the main steps to extract the 
divergent part. We start with the integral representation: 



In = 



1 



ds 



du ■ 



(34) 



where = ]\P + 2fe3U'^/A'*, Aw^ = 2biu + 262U^/A^. Expanding the exponential in powers of Aw^, it is easy to 
see that only the leading term is divergent for Z) — > 3 and can be explicitly evaluated. Then, performing a series 
expansion on e, we can write 



A2 



- + ln(M) 



finite + 0{e), 



where the finite terms arc independent of both /i and e. Inserting this in Eq. (j33p . we choose the counterterms as 

1 



(35) 



5m 



5\2p = - 



so to cancel only the poles in e = 0. Afterwards, the e 
equation reads 



87r2 

^2p+2 



87r2 



2&3 e 
A2 1 



(36) 



(37) 



limit can be safely taken, and the renormalized evolution 



□l(x) 



AfA^ 
87r2V26; 



In (/i) — finite j 00 — 



\R A 2 




In (/i) + finite 



(2p-l) 



= 0. 



(38) 



This equation inherits a ^ dependence from the renormalization (or regularization) procedure, which remains even 
after taking the e — >■ limit. However, the physics must be independent of the scale so introduced. Therefore, the 
ln(/x) derivative of the previous equation must vanish term- by-term in the sum of powers of 0o • Then, it follows that 
renormalized couplings must satisfy the following UV RG equations: 



dX2p 

dfj. 



AfA2 



87r2 

^2p+2 



2&3 

A2 



87r2 



/2fe3 

ll2\ 



(39) 
(40) 



To conclude this section, we introduce the notation {(j)^)ren = {4'^) — for the finite part of (02^'4_d(o) 

the contributions of higher adiabatic order. Then, the renormalized evolution equation is 




0. 



(41) 
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This expression will be picked up later in Sec. I VII to discuss the type of couplings with the curvature that come from 
the second-adiabiatic-order part of {(jp)ren- 

V. INFRARED REGIME 



The minimal subtraction scheme cannot be used to study the infrared behavior of the running couplings, but a 
more physical renormalization scheme must be used instead, such as the momentum subtraction scheme. Another 
possibility to study this problem is to use the so-called RG-flow equation for the effective action [l^. As we will see, 
this approach is useful not only for obtaining the RG equations for the coupling constants of the potential in the UV 
and IR regimes, but also to extend the analysis beyond the one-loop approximation. 

We will therefore generalize here the methods presented in Ref. [13] to the Lifshitz case. In this formalism, the 
1-loop correction to the effective potential can be written as 



2n 
1 

2 



d(l){x)d(j){y) 



(42) 




fc2 + y"($) -\TL{k^) 



(43) 



where the integral is performed over the Euclidean momentum space k. The renormalization scale ^ is identified with 
the IR momentum cutoff, so that the integration region is the volume contained between two concentric 3-spheres of 
radii and A[/y. 

The main difference with the Lifshitz case is the anisotropic scaling. This is evidenced in the integrand, which, 
with the modified dispertion relation (|24l) . gets the replacement 



-kl + Co''{k) = kl+u^{k), 



(44) 



with uj'^{k) ~ ijj'^{k) — A'P [the effective mass term is included in the potential X^($)], and then it is no longer 
spherically symmetric in 4 dimensions. Hence, the anisotropic scaling should be taken into account by modifying also 
the integration region, using different cutoffs in the spatial and temporal directions. 

We motivate one way of performing this generalization by first looking at the integral over the 3-sphere from the 
usual Lorentz-invariant case. This integral can be written 

/ dkk^= dkk^- dkk^= dk \k\^ dk4 

Jfi Jo Jo Jo J -J J^lv-\^? 



dk / dk4. 

J-y/iJ.^-\k\^ 



(45) 



In this form, we can introduce the anisotropy by modifying the limits of the k^ integral in the following manner: 



A2 



\k\' 



CJ■^{^l)-Co■^{k). 



(46) 



This leads to an integral over the volume contained between two sort of prolate 3-spheroids, with their major axes 
aligned in the k^ direction, and whose cross sections in the plane perpendicular to k^ are concentric 2-spheres. This 
3- "rugby-ball" can be continuously deformed into a 3-sphere when w^(fc) — >■ |fc|^, which is expected to happen in the 
IR. 

Putting all this together, and taking into account that the integrand now only depends on |fc| and the expression 
is generalized to 



^ ^ ' 4^3 .;o 



47r3 



dk \k[ 



dk IfcP 



V"^(At;v)-<i2(fc) 



dk^ In 



V (<&) 

u^k) 



dkd In 



kl+Q-2{k) 



(47) 
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We cannot go much further from here to obtain JJ^"'^^ ($) because of calculational difficulties. Nevertheless, in order to 
calculate the RG equations, we can take the ln(/i) derivative of uli'\^) before performing the integrals, which greatly 
simplifies the computations. 

The renormalized constants arc obtained by taking a given number of <f> derivatives of C/^(^') and then evaluating 
at $ = 0: 



and, in general, 



(48) 
(49) 

(50) 



From the expression (j47p . the first term is independent of fi and is dropped, while in the second, there are two nested 
integrals whose limits depend on /i. Performing the derivatives, we get 



dii 



/i du 
47r3 dfc^ 



dk IfcP In 



(51) 



On the right-hand side, the logarithm becomes independent of |fc| once we evaluate 
can take it out of the integral. What remains is quite simple: 



ijo'^in) —u]'^{k). Hence, we 



M = -f{^^dh}) In 



1 



(52) 



where we have defined 



/(M,W) = ^^^ / dk 



Co^{li)-Cj^{k) 



(53) 



For the dispertion relation ui'^{k) ~ 2bik^ + 2&2fc"'/A^ + 2&3fc^/A'*, this function /(/^, {&i}) cannot be calculated 
analytically, making it difficult to reach expressions valid for any value of fi. However, it is possible to study its 
asymptotic behavior in the UV and IR limits. On the one hand, in the UV we have A and therefore uP'{k) ~ ^k^, 
which leads to 



fuv - fip-, h) 



87r2 A2' 

while on the other hand, in the IR, /i ^ A so that uj'^{k) ~ 2bik^, giving 



(54) 



fiR - /(m,&i) 



26i 4 



(55) 



Up to here, we have kept the 1-loop corrections to the effective potential. A RG improvement [o] can be made on 
Eq. (|5^ by replacing V (<f>) by 9|,C/^ on the right-hand side, leading to 



M = -./(a*, V>,}) In 



1 



(56) 



When taking the $ derivatives and evaluating at $ = to find the RG equations for m|j,, A^^ or any A^, only the 
logarithm in the previous expression comes into play. Within the resulting expressions, there will be terms containing 
different number of $ derivatives of f/^, which might vanish upon evaluation at $ = 0. For instance, for the potential 
considered here (jlOp . odd derivatives of C/^ will evaluate to zero, as also will any number of derivatives greater than 
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We start by considering the RG equation for the renormalized mass m|j. 
evaluation at $ = of Eq. ([55)1 , it reads 



After the appropriate derivation and 



-fi^^Ab^}) 



A4^ 



(57) 



The asymptotic behavior of this equation is obtained by replacing tD^(/i) and /(/i, by their corresponding ex- 
pressions in the UV or in the IR. The results are 



dfi 
dm^ 
dfi 



AfA2 



167r2 



A-^fi (UV), 
Ai<A (IR). 



(58) 
(59) 



While the first equation is the same that was obtained earlier by dimensional regularization and minimal substraction 
([55)1 . the second equation reproduces the standard result for the Lorentz- invariant case for bi — >■ 1/2. 

In a similar fashion, we obtain the RG equation for A^^ by taking four derivatives prior to the evaluation at $ = 0. 
The result is 



M^ = -/(/^,{6.}) 
dn 



3(Af)2 



(m| + w2(^)) (r74 + c:;2(^)) 



(60) 



and the corresponding UV and IR asymptotic expressions are 



dXf 
d/i 

dXl 
djjL 



AfA^ 



3(Af)2 



A 

8^2^25^ ' 8^2(253)3/2 y- 
167r2 {ml + 2biy?) 



(Af)V^ 



167r2 



2bl^l^Y 



A</i (UV), 
At<A (IR). 



(61) 
(62) 



As the second term in the right-hand side of the UV equation can be neglected in this limit, it again coincides with 
the one found previously pO)) for p = 2, while for the IR equation, we recover the usual case provided bi — )- 1/2. 

So far, we have dealt with the mass and a coupling that are both renormalizable in the usual case. We would like 
to see what happens with the new couplings that are not renormalizable in that case. However, for couplings with 
p > 2, the corresponding expressions get more and more complicated as a higher number of $ derivatives are taken, 
generating increasingly more terms. This can be thought diagramatically as the increasing number of ways to obtain 
a 2p-legged diagram at 1 loop for higher p's, because there are more types of vertices with fewer legs available. 

Nevertheless, let us consider Ag^ as the easiest example of a coupling that it is not renormalizable in the usual case. 
Its RG equation is 



dA£ 



-/(m,{&J) 



(^2(^)) 



15AfA« 



(m| 



■'^^(m)) 



30(Af)3 



(63) 



and the UV and IR limits are, respectively. 



dA£ 
dfi 

dxl 

dfi 



AfA2 



15A|^A 



R\R 
6 



A 



87r2y25:^ 87r2(263)3/2 Vm 
V2h Al^^i* 157261 



30(AiO 



R\3 



8^2(263)5/2^2 

A4 Ag [I 



167r2 (m 
30^/26T 



+ 26i/i2) 

(A4«)3^4 



167r2 



167r2 



{ml 



2bi^i^Y 
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A<^i (UV), 



At<A (IR). 



(64) 



(65) 



Once more in the UV equation, the second and third terms in the right-hand side can be dropped in this limit, and 
we get the same as with dimensional regularization and minimal substraction, that is Ea. (|40p for p = 3. In the IR 
limit, the equation is the usually encountered in effective field theories at low energies, again for bi — > 1/2. 

Higher order couplings get more and more complicated expressions. A full analysis of the running of the coupling 
constants would involve a solution of the coupled differential equations for all couplings. 
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VI. COUPLINGS TO THE CURVATURE 



We will now consider the coupling of the scalar field to terms containing second derivatives of the metric. In the 
z = 1 case, {(f)'^) contains a divergence proportional to the Ricci scalar, and therefore the corresponding counterterm 
^'^'^'>R(jP should be introduced in the action. For z = 3, the second adiabatic order correction ([30|) is finite. These 
terms can be grouped to form the linearized expressions for the following curvature invariants: 



K 



h - 2diSN'' 



(66) 
(67) 



where ^^^R is the scalar of 3-curvature, and K is the trace of the extrinsic curvature Kij. The dot stands for time 
derivative. With these identifications and using the following relation valid at linear order on the metric perturbations 



(4) 



12.AD{2) ^ A''^^R + Bk + c{-d^5n-2k 



we can rewrite Eq. ([50)) as 



Here, the coefficients A, B and C come from the integrals defined in Eq. ([51]) . and their expressions are 

1 r . a;-i/2 



A 



B 



C 



967r2 



dx 



(1 - 26i 
487r2 

1 



„l/2 



dx 



967r2 



dx 







[Af2/A2 + 2bix + 2&2x2 + 263a;3]3/2 ' 

x3/2(452 +663X) 

[M2/A2 + 2bix + 2b2X^ + 2b3X^f^'^ ' 



(68) 
(69) 

(70) 
(71) 
(72) 



The renormalized evolution equation (|4ip then reads 



2 



A + Bk + C[ -d'^5n - 2k 




{2p-l)\ 



(73) 



Had wc taken the parameters bi = 1/2, &2 = &3 = 0, or A — >■ 00 from the beginning, the coefficients B and C would 
be identically zero, and the correction (j69p would have been equal to the standard result for the Lorentz-invariant 
theory (see Ref. [22|), that is 



(74) 



where Ausuai is A with its integrand evaluated for these values of the parameters before integration, and it is a 
divergent integral. However, if wc consider a Lorentz-noninvariant UV completion of the theory like the one wc are 
studying in this work, the situation is different. When taking the infrared limit, which is A — >■ 00, the coefficients A, 
B and C do not vanish, not even in the particular case 61 = 1/2, 62, 63 — )■ 0. 

On the one hand, the coefficient B can take various values depending on &i's initial value. If bi is precisely 1/2, 
then _B = 0, but for any other case, it will depend on the particular running of bi, 62 and 63, with _B 00 being a 
possibility. To study this, the running of bi should be studied in more detail. 

On the other hand, the coefficient C is given by an integral that is finite regardless of the mass M, so it can be 
calculated explicitly by taking the A — >■ 00 limit (or M — > 0) prior to integration. The result is finite. 



C = 



V2{2b2 + iVhh) 
Vh{b2 + 2Vhh) ' 



(75) 
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and has nonvanishing limits even when 62 — >■ and 63 —> 0, effectively leaving a non generally-covariant remnant in 
the 1-loop evolution equation of (0) that does not vanish in the IR, in spite of the parameters taking their Lorentz- 
invariant values. Similar effects have been observed before, for example, in the case of a theory with a Yukawa-type 
scalar-fermion interaction in flat spacetime and with a generically modified dispersion relation for the fermion at high 
energies (23j . 

We will now take a closer look to the integral (fTO]). that defines the A coefficient to try to find out its behavior 
when A — 00. In this limit, the integral has an infrared divergence; however, away from it, k = Af/A acts as an 
infrared cutoff. By noting that the integrand is sensitive to the value of k = M /A only when a; <C 1 , we can separate 
the integration interval in two parts, < x < e <C 1 and e < x < 00. In the first region, we can drop and a;"* in 
front of x'^, while in the second region, we no longer have the infrared divergence and, consequently, we can safely 
take the limit k — > 0. Then the coefficient A approximately reads 



1 



967r2 



dx 



1 



dx ■ 



which now contains two integrals that can be performed analytically, giving 



A 



1 



967r2 



1 
71 



ArcCosh 1 



45ie 



In 



2bi + bze + 2v/bi(&i + hit + 636^) 
626 + 2\/M^e 



(76) 



(77) 



The dependence in k, and ultimately on A, is contained in the first part. As we are interested in the behavior for 
K — >■ 0, we perform a series expansion in k around zero. Meanwhile, the second part is finite and of order 1 for any 
positive sufficiently small nonvanishing value of e, which is when this treatment gives a closer estimate of A. Finally, 
we can say that the behavior of this coefficient for large A is approximately 



A: 



1 



967r2VFr 



V21n 



0{\) 



(78) 



which shows that the constant A, divergent in the Lorentz-invariant theory, is now finite with A acting as a cutoff. 
The dependence is logarithmic, which means that the corrections are not large unless M is many orders of magnitude 
below A. 



A. Coupling to "healthy" terms 



According to Ref. |18j . the Hofava-Lifshitz theory of gravity is better behaved when invariants formed with 
ai = di\a{N^ are included in the gravitational action than when they are not. The most relevant invariant is a^a*, 
which is absent in general relativity. 

It is interesting to see whether the self-interaction of the scalar field induces couplings to such terms or not. In 
order to do that, it is necessary to go beyond the weak field approximation considered so far. Instead of tackling this 
complicated calculation, we will follow a shortcut by performing a derivative expansion on A^, taking as a starting 
point a metric of the form 

ds^ ^ ~N^{x)dt^ + dx^ (79) 



(80) 



which is enough for our present purpose. The corresponding equation for the Green function is 

-^d^Gix, x') - NM^G{x, x') + 2&i9, {Nd'Gix, x')) 

-^d^ {Nd^Gix, x')) + {Nd^Gix, x')) + {Nd^G{x, x')) = d{x - x'). 

After expanding in spatial derivatives of N up to second order and then going to Fourier space, the equation reads 

dG{k) 



iV2 



G(fc) 



ftp 



1 

+2 



iV2 



2\ d'Gik) 



d{k'^Y 



k^5'^G{k) 



d{BY 



-k'k^ 




aG(fc)\ 

dh J 



dk,dk. 



(81) 
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To solve it, the solution is written as a sum of contributions of ascending order in derivatives of N, noted with a 
subscript Ni, so 



G(fc) = Gjvo(fc) + G^i(fc) + Gjv2(fc) 



(82) 



and then the equation is solved iteratively order by order. The expressions obtained for GNQ^k), Gj\ji{k) and Gj\j2ik) 
are given in the Appendix. 

The calculation continues from this point in the same manner as in the weak gravitational field approximation, 
giving the following results: 



)N0 



)N1 



/N2 



87r2' 
0, 

1 



967r2 



2/o + 

6 



1^' 



(83) 
(84) 

(85) 



where the integrals /q, Iq, I2 and I3 are the same as defined earlier in Eq. pip . Notably, the coefiicient accompanying 
a^a* has vanished, so the result is just the same as the one obtained in the previous section (pO| . 

Evidently, the absence of this kind of correction is more fundamental than a simple over-approximation. In fact, it 
can be easily verified that the result would be the same if we had taken the Lorentz-breaking terms in the action to 
be slightly different, as we will now see. 



B. More general Lagrangians 



We already found in the previous section that couplings between the scalar field and 0^0* seem not to be generated. 
However, one can ask if this is a consequence of the particular form of the action we started with Eq. (fTO|) . Indeed, 
we could have chosen any term of the form 



(86) 



where the superscripts a and b denote the number of spatial derivatives acting on the field which satisfy a+b = 2n < z, 
regardless of how the indexes of each pair of derivatives are contracted. Performing a variation of the action, this 
kind of term leads to 



5 



(fxdtNd''(bd^ 



d^xdt [(-l)"9"(iVa^0) + {-lfd\Nd''(j))\ 
= {-If [ d^xdt \2Nd''+''(t> + (a + b)d,Nd'd''+''-\ 



L[AN, ( 



(-1)" / d^xdt 



2Nd^''6 + 2nd,Nd'd 



i f)2n — l. 



+ L[AN,< 



(87) 



where L contains terms with second-order-or-higher spatial derivatives of N, and its form will be dependent on a 
and b. Nonetheless, the two terms relevant for this discussion do not dependent on a nor b, but on the sum of both. 
Therefore, it is irrelevant how we choose to distribute the derivatives; the result will continue to be the same, i.e. 
terms proportional to 0^0* will not be generated. 

We can also argue how the results will change if we consider interactions with derivatives of the field in the 
action. For example, if we add a term of the type N^cj)^ Acj), the 1-loop evolution equation of (0) will get corrections 
proportional to {(j>A(j)), {(f>dicf>) and {dicfid'^ (j>) , besides the already known (0^). These are calculated from the imaginary 
part of the Feynman propagator, taking the appropriate derivatives prior to taking the coincidence limit. This will 
bring down powers of fc in a Fourier space representation, so the second adiabatic order part of these expectation 
values can be expressed as 



d'^k 



fik')GN2ik), 



where Gjv2(fc) is the same from Eq. ((5^ and is given in the Appendix. As it happens in the calculation of Eq. 
the terms proportional to a^a* vanish before the d'^k integrals are performed. Hence, the presence of these extra fc' 
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factors does not change that fact. Accordingly, interaction terms with derivatives of the field seem not to generate 
corrections with a^a' either. 

It remains to be shown whether it is true in general that Lagrangians of the form 

^ = N^C, (89) 

with C 7^ C{diN), do not generate quantum corrections that couple {(j)) with a^a*. On the other hand, these corrections 
are generated if the classical Lagrangian contains explicit couplings of the quantum fields with diN, see for instance 
Ref. HI. 

VII. CONCLUSIONS 

In this paper, we studied some aspects of Lifshitz-type field theories in curved spaces, towards the main goal of 
assessing the phenomenological viability of theories with broken Lorentz invariance. We focused on two aspects: 1) 
the development of methods for analyzing the running of the constants in the Lagrangian (and with that, the eventual 
recovery of the Lorentz invariance at low energies); 2) the study of nonminimal couplings between matter fields and the 
spacetime metric. Regarding the second point, we looked at two types of terms: couplings with the spatial derivative 
of the lapse function (i.e., with 0^0' where = dilnN) and covariant and noncovariant couplings with the curvature 
scalars. 

We presented a detailed study of the 1-loop renormalizability of z = 3 self-interacting Lifshitz fields in curved 
spacetimes, for interactions of the form As the higher spatial derivatives of the propagators improve the ultraviolet 
behavior of the theory, the divergences appear only in the zeroth adiabatic order; that is, they are essentially the flat 
spacetime divergences. For the computation of the RG equations for the coupling constants both in the UV and in the 
IR, we generalized the exact RG methods to the Lifshitz case. A full analysis of the running of these couplings would 
require us to solve a complicated system of coupled nonlinear differential equations. However, a more interesting point 
would be to analyze the running of the constants corresponding to the kinetic terms. Therefore, it would be worth 
it to extend the method of the exact RG we have implemented for the effective potential to the effective action in 
order to study the IR running of derivative couplings [16|. For this, one could use a generalization of the derivative 
expansion methods described in Refs. (25l [26} . This would be a very useful tool for analyzing the running of the 
couplings corresponding to the IR-relevant terms with two (time and spatial) derivatives and to study the eventual 
emergence of a universal limiting speed. 

We also computed the second adiabatic order contributions to the mean value equation of the Lifshitz field to 
analyze the couplings to the curvature. We found that the self-interaction induces couplings to the four-dimensional 
Ricci scalar, the 3-curvature and to the extrinsic curvature, but it does not generate couplings with a^a'. We argued 
that the absence of quantum corrections involving such derivatives of the lapse function seems to be general, provided 
the Lagrangian is of the form given in Eq. (j89|) . 

The coupling with the four-dimensional Ricci scalar diverges in the limit A —> 00, as expected from the fact that 
it diverges in the usual theory. However, the divergence is logarithmic in M/A. Therefore, unless the effective mass 
M is many orders of magnitude below A, in these theories, the coupling to the four-dimensional Ricci scalar does not 
receive large quantum corrections ^ 3> 1. The noncovariant couplings to the 3-curvature and the extrinsic curvature 
are expected since general covariance is broken already at the classical level. However, it is remarkable that there is 
a nonvanishing finite remnant even in the limit A — > cx3. 

The nonminimal couplings to the curvature may have interesting consequences both in the early Universe and in 
astrophysics, in regions of strong gravitational fields. For example, there are proposals in which the Higgs boson is 
supposed to be responsible for infiation [27j . In the standard (z = 1) scenario, the key point of this proposal is the 
nonminimal coupling of the Higgs field to gravity It has been shown that ^ 10"* is required in order to suc- 
cessfully obtain infiation and a power spectrum of the primordial perturbations in good agreement with observational 
data [Higgs-infiationary models with realistic quadratic Higgs self-interaction and minimal coupling (^ = 0) generates 
an unacceptably large power spectrum of the primordial perturbations] [l^l- From our results, it is clear that such 
larg e couplings could be generated if M <|; A, or in the presence of a mechanism similar to that proposed in Ref. 
|15| . in order to enhance the logarithmic dependence. Like the usual coupling with ^'^^R, the noncovariant ones could 
also be relevant in inflationary models where a Lifshitz field generates the primordial perturbations 0, • These 
proposals have been questioned on the ground that they might suffer from unitary problems, due to the large value 
of the coupling to the curvature (see, for instance, Refs. [29l - l33j ). Recently, in the context of Hofava-Lifshitz gravity, 
it has been shown that infiationary models with a z = 3 "Higgs" field with ^ = could be constructed [l^, so that 
one would not need to worry about preventing unitarity problems. It would be interesting to further analyze these 
kinds of models. 
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For a nonminimally coupled z = 1 Higgs field, observational consequences in strong-gravity astrophysical environ- 
ments of the curvature dependence of the mass generation have been analyzed in Refs. [SJ, H^- As ^''^i? vanishes 
for spherically symmetric sources described by the Schwarzschild metric, in order to obtain nonvanishing effects, the 
author assumed an unnatural coupling to the square root of the Kretschmann invariant Ki = R^^p^rR^'^^'^ . For 
the case of z = 3 Lifshitz fields, the interactions generate couplings to ^'^•'i? and K, which do not vanish for the 
Schwarzschild metric. Therefore, it is not necessary to assume a coupling to the Kretschmann invariant. In Refs. 
[U,!!^ it was found that, unless one allows quite large values for the associated coupling constant to y/Ki, the effects 
are quite far from what can be achieved with any foreseeable survey. Hence, we can expect a similar conclusion for 
the noncovariant couplings. 

Finally, as for future work, it would be interesting to consider theories containing fermions and gauge fields. From 
the phenomenological point of view, in addition to studies of the running of the coupling constants in flat space, it 
would also be of interest to analyze the nonminimal couplings induced by quantum corrections — in particular in the 
context of electrodynamics in curved spaces — and use astrophysical and cosmological observational data to constraint 
the parameters of the action, generalizing the analysis of Ref. [s^ to the case of noncovariant couplings to the 
spacetime metric. 
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Appendix A: First and second adiabatic order expressions for the propagator 

In this Appendix, we provide some explicit expressions for the adiabatic expansion of the propagator. In the 
case of the weak gravitational field described in Sec. HI, the first order in this approximation involves the following 
expression: 



fk{p) = [i^-Sh]{kl + kopo) + {poh + 2kok, + piko)6N' - [-+Sfi]ul 
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(90) 
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Then, the first and second adiabatic orders for the propagator are given by 



- Dl{k)k^diy^kjkk 
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(92) 



On the other hand, for the metric considered in Sec. VI A, the expansion of the propagator in derivatives of the 
lapse function N is 
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where fco = ko/N. 
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